ONR  Graduate  Traineeship  Award  in  Ocean  Acoustics  for  Sunwoong  Lee 


PI:  Prof.  Nicholas  C.  Makris 
Massachusetts  Institute  of  Technology 
77  Massachusetts  Avenue,  Room  5-212 
Cambridge,  MA  02139 

phone:  (617)  258-6104  fax:(617)253-2350  email:  makris@mit.edu 

Co-PI:  Sunwoong  Lee 
Massachusetts  Institute  of  Technology 
77  Massachusetts  Avenue,  Room  5-435 
Cambridge,  MA  02139 

phone:  (617)  253-1032  fax:  (617)  253-2350  email:  sunwoong@mit.edu 
Grant  Number:  N0001405 10073 


LONG-TERM  GOALS 

The  long-term  goal  of  this  research  is  to  develop  optimal  array  signal  processing  techniques  for 
passive  source  localization  in  littoral  shallow-water  environments.  It  has  long  been  known  that  multi¬ 
modal  dispersion  in  a  shallow  water  waveguide  degrades  the  perfonnance  of  bearing  estimates  by 
conventional  plane-wave  beamfonning.  This  is  due  to  spurious  effects  unique  to  the  waveguide 
environment  such  as  multiple  peaks  and  beam  spreading  in  the  beamformer  output  [1,2].  We  have 
developed  array  signal  processing  techniques  that  account  for  and  exploit  multi-modal  wave 
propagation  and  dispersion. 

OBJECTIVES 

The  primary  objectives  of  the  research  are  to: 

•  Develop  the  array  invariant  passive  localization  techniques  that  require  little  a  priori 
knowledge  of  the  wave  propagation  environment. 

•  Instantaneously  and  simultaneously  localize  multiple  sources  in  a  shallow-water  waveguide 
without  extensive  computations  and  without  ambiguity. 

APPROACH 

The  array  invariant  method  for  passive  localization  of  impulsive  sources  has  been  theoretically  derived 
and  experimentally  demonstrated  in  FY05  effort  using  the  data  acquired  form  the  Acoustic  Clutter 
Experiment  in  2001  and  2003.  It  has  been  shown  that  simple  and  robust  passive  source  localization  is 
possible  using  the  array  invariant  method  without  a  priori  knowledge  of  environmental  parameters. 

We  have  extended  the  array  invariant  method  for  range  and  bearing  estimation  of  a  broadband  random 
noise  source  that  is  not  necessarily  impulsive  in  the  time  domain,  while  maintaining  the  advantages  of 
the  array  invariant  method  for  impulsive  source  localization.  We  have  theoretically  and  numerically 
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shown  that  the  array  invariant  method  can  be  used  not  only  for  localization  of  a  single  source,  but  for 
simultaneous  localization  of  multiple  noise  sources  in  an  ocean  waveguide  without  ambiguity. 

WORK  COMPLETED 

The  research  work  in  Fiscal  Year  2006  was  a  great  success  in  extending  the  array  invariant  source 
localization  techniques  in  a  dispersive  waveguide  that  require  little  a  priori  knowledge  of  the 
environment.  The  extended  array  invariant  techniques  enable  instantaneous  and  simultaneous 
localization  of  multiple  random  noise  sources  in  a  horizontally-stratified  ocean  waveguide  from 
passive  beam-time  intensity  data  obtained  after  conventional  plane-wave  beamforming  of  acoustic 
array  measurements. 

Localization  of  multiple  sources  using  the  array  invariant  method  has  advantages  over  the 
conventional  triangulation  method  using  two  arrays.  First,  the  array  invariant  method  does  not  require 
two  arrays  with  sufficient  spatial  separation,  which  is  a  significant  advantage  when  the  localization  is 
made  from  a  mobile  platform.  Second,  the  array  invariant  method  does  not  suffer  from  the  ghost 
source  problem  typical  in  the  triangulation  method  when  there  are  multiple  sources  present  [3].  The 
array  invariant  method  is  significantly  advantageous  over  MFP,  where  unambiguous  source  range 
estimation  is  nearly  impossible  even  with  accurate  environmental  knowledge,  when  there  are  unknown 
number  of  multiple,  spatially  stationary  sources  in  an  ocean  waveguide  [4,  5] . 

RESULTS 

We  have  previously  shown  that  the  localization  of  an  impulsive  source  can  be  achieved  using  the  array 
invariant  method  without  requiring  the  a  priori  knowledge  of  environmental  parameters  and  without 
extensive  computations  [6].  We  have  further  extended  the  method  to  show  that  the  array  invariant 
method  can  be  applied  for  localization  of  a  broadband  random  noise  sources.  This  is  achieved  by 
cross-correlating  the  beam-time  series  with  the  time  series  measured  by  the  acoustic  sensor  at  the 
center  of  the  array  [7]. 

We  first  show  that  the  range  and  bearing  of  a  broadband  random  noise  source  can  be  determined  using 
the  extended  array  invariant  method.  The  source  is  located  at  rQ  =  5  km  and  0o  =  60°  as  shown  in  Fig. 

1.  The  environment  is  100-m  deep  Pekeris  waveguide  with  sand  bottom.  The  source  and  the  receiving 
array  depths  are  50  m  and  30  m,  respectively.  The  aperture  of  the  receiving  array  is  150  m.  The  source 
power  spectral  density  is  assumed  to  be  0  dB  re  1  pPa  /  Hz  within  the  390  to  440-Hz  frequency  band. 
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Figure  1:  Top  view  of  the  geometry  for  the  single  source  localization  example.  The 
source  is  at  60°  from  the  broadside  of  the  horizontal  receiving  array, 

and  at  5-km  range. 
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Figure  2:  The  cross-correlated  intensity  field  for  the  source-receiver  geometry  shown 
in  Fig.  1  in  the  sand-bottom  Pekeris  waveguide  environment.  The  black  solid  line 
overlain  in  the  t  >  0  domain  is  the  beamformer  migration  line  in  terms  of  the  reduced 
travel  time.  The  black  solid  line  in  the  r  <0  domain  is  at  sin  90 . 


The  cross-correlated  intensity  field  IBo(s,t)  is  shown  in  Fig.  2,  where  s  -  sin 6 is  the  array  scan  angle, 

and  r  is  the  delay  time.  Two  black  solid  lines  are  overlain  on  Fig.  2.  One  line  in  the  r  >  0  domain  is  the 
beamformer  migration  line  ?(r)  for  an  impulsive  source  in  an  ideal  waveguide  as  defined  in  Eq.  (8)  of 


Ref.  [6],  and  the  other  line  in  the  r  <  0  domain  marks  the  source  bearing  sin  6o .  These  two  lines  show 

that  the  cross  correlated  intensity  field  is  bounded  by  these  lines  except  that  the  field  is  smeared  by  the 
beampattern  of  the  horizontal  receiving  line  array.  Given  these  two  bounding  lines  in  the  cross- 
correlated  intensity  field,  we  can  estimate  both  source  range  and  bearing  without  ambiguity. 

In  order  to  identify  these  two  bounding  lines  from  Fig.  2,  we  apply  an  image  transfonn  technique 
similar  to  the  Radon  transform  [8]  to  the  cross-correlated  beam-time  image.  In  this  transfonn,  a  given 
image  is  integrated  along  a  semi-infinite  line  starting  from  (5,  0)  and  having  an  angle  <j>r  from  the 
positive  x  -axis  in  the  clockwise  direction.  Now  we  show  that  both  the  range  and  bearing  of  the  source 
can  be  extracted  from  the  transformed  intensity  field  in  Fig.  3. 


s  =  sin  0 

Figure  3:  The  transformed  intensity  image  In(s,<f>r)  of  the  cross-correlated  intensity 

image  in  Fig.  2. 

First,  the  bearing  of  the  source  can  be  estimated  from  the  value  of  $s$  where  the  maximum  of 
I„(s,</>r  =  180°)  occurs.  In  Fig.  3,  this  peak  is  seen  to  occur  at  sa  =  0.853  or  0o  =  58.5°  .  This  is  within 
3\%  of  the  true  source  bearing  0o  =  60° .  Second,  once  the  bearing  of  the  source  has  been  estimated, 
the  range  of  the  source  can  be  estimated  using 

„  _  c(z)  sin  #o 
tan$. 

where  <j>r  =  argmax^  Ifs  =  s,<f>r ).  In  Fig.  3,  this  peak  occurs  at  <f)r  =  346.7°  .  The  range  of  the  source  is 
then  estimated  to  be  ra  =5.4  km,  which  is  within  10%  of  the  true  source  range. 


The  approach  used  for  localization  of  a  single  source  can  also  be  used  for  localization  of  multiple 
uncorrelated  random  noise  sources  in  an  ocean  waveguide  in  both  range  and  bearing.  Here  we 
consider  the  case  where  there  are  3  uncorrelated  noise  sources  Si,  S2,  S3,  as  shown  in  Fig.  4.  The 
source  Si  is  at  range  r/  =  8.2  km  and  bearing  $\theta_l  =  45A\circ$.  The  other  two  sources  S2  and  S3 
are  at  the  same  bearing  $\theta_2  =  \theta_3  =  60A\circ$,  and  their  ranges  are  r2  =  3  km  and  r3  =  10  km, 
respectively.  We  assume  that  these  three  sources  have  the  same  source  power  spectrum  of  0  dB  re  1 
|iPa2  /  Hz  in  the  390  to  440  Hz  frequency  band,  but  they  are  uncorrelated  with  each  other. 
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Figure  4:  Top  view  of  the  geometry  for  the  multiple  source  localization  example.  All  3 
sources  are  at  30-m  depth,  and  the  receiver  array  depth  is  50  m.  The  receiver  array 

aperture  is  150  m. 


The  cross-correlated  intensity  field  IBo(s,t)  for  this  multiple  source  scenario  is  shown  in  Fig.  5.  The 

cross-correlated  intensity  in  Fig.  5  exhibits  two  distinct  intensity  fields  when  r  <  0,  one  at  the  bearing 
of  Si,  and  the  other  at  the  bearing  of  58  and  S3.  When  r  >  0,  the  field  for  55  and  S3  splits  into  two 
distinct  fields,  since  their  ranges  from  the  receiver  array  are  different.  The  cross-correlated  intensity 
field  for  Si  when  x  >  0  shows  the  same  pattern  as  that  of  the  single  source  scenario. 
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Figure  5:  The  cross-correlated  intensity  field  l  Ro(s,z)  for  the  source-receiver  geometry 
shown  in  Fig.  4.  The  two  black  solid  lines  overlain  in  the  r  <0  domain  marks 
sin  0]  and  sin  02 ,  respectively.  The  three  black  solid  lines  in  the  r  >  0  domain  are  the 
beamformer  migration  lines  for  sources  Si,  S2,  and  S3  in  terms  of  the 

reduced  travel  time. 

We  now  can  apply  the  image  transfonn  method  previously  introduced  for  localization  of  single  source, 
and  localize  all  three  sources  in  range  and  bearing  without  ambiguity.  The  detailed  procedures  and 
estimation  results  are  provided  in  Ref.  [7]. 

IMPACT/APPLICATIONS 

The  array  invariant  method  developed  from  Fiscal  Year  2006  effort  has  significant  advantages  over 
existing  source  localization  methods  for  practical  source  localization  scenarios.  This  is  because  the 
array  invariant  method  does  not  require  a  priori  knowledge  of  the  environmental  parameters,  nor  does 
it  require  extensive  computations. 

The  array  invariant  method  has  tremendous  potential  application  in  shallow-water  surveillance 
missions  and  anti-submarine  warfare.  The  method  enables  instantaneous  and  simultaneous  localization 
of  multiple  broadband  noise  sources  in  a  littoral  environment  using  towed  arrays  from  surface  ships  or 
submarines. 
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